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a) Objectives: 

The task of algorithm-development activities at USF continues. The algorithm for 
determining chlorophyll a concentration, [Chi a ] and gelbstoff absorption coefficient for 
SeaWiFS and MODIS-N radiance data is our current priority. 


b) Task Accomplished: 

A preliminary algorithm for chlorophyll a and gelbstoff parameterized for the Gulf of 
Mexico in summer has been accomplished. The bio-optical algorithm being developed for 
SeaWiFS and MODIS-N is based on a semi-analytical model of remote-sensing reflectance 
(R rs ) and is an extension of the algorithm discussed in Carder et al. (1991). The major 
difference between the two occurs in the spectral term for the effects of backscattering and 
the upwelling distribution of radiance as developed in Lee et al. (1994). As a result of this 
change a bio-optical algorithm utilizing the SeaWiFS wavebands centered at 412, 443, 490, 
555, and 670 nm is used to estimate chlorophyll a concentration ([Chi a ]) and gelbstoff 
concentration, represented by its absorption coefficient at 400 nm (a g (400)). 

The R rs model has numerous parameters that cannot be fixed and applied to the entire 
globe; i.e., they are site- and season-specific. For example, absorption per unit chlorophyll 
by phytoplankton can change with species, and with nutrient and lighting conditions by as 
much as a factor of five (Morel and Bricaud 1981; Carder et al. 1991; Morel et al. 1993). 

Also, particle size and concentration both have a significant effect on the spectral 
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backscattering coefficient, b b (X), of ocean water, because pure water backscatters as ~ X" 4 , 
large particles backscatter as ~ X* 0 , and smaller diameter detritus and bacteria backscatter 
with a spectral dependence between the two (Morel and Ahn 1990, 1991). If many of these 
factors did not covary, the simple wavelength-ratio algorithms of the CZCS (Gordon and 
Morel 1983) would never have worked as well as they did. In trying to understand these 
covariances, we have developed semi-empirical expressions that depend on chlorophyll for 
some components of the model. 

Extensive field data sets are needed to allow seamless modification of the model 
parameters with time and space. The changes required will be due mostly to changes in the 
dominant plankton groups present and the subsequent effects on bio-optical parameters such 
as pigment packaging. Acquiring such data sets on a global scale should be a major 
community goal during the next few years. We have developed a scenario that can both 
guide the parameterization process and provide an initial regional implementation of the 
algorithm. 

Two cruises to the Gulf of Mexico were conducted, one during the spring and one 
during the summer of 1993, to acquire data for development and testing of the algorithm. 
Pigment data from the spring cruise are not yet available for model validation, but the optical 
and bio-optical portions of the algorithm have both been tested on the summer data, where 
[Chi a] ranged from < 0.10 to > 40 mg m' 3 , and gelbstoff-to-phytoplankton absorption 
ratios at 443 nm (a g (443)/a^(443)) ranged from 0.3 to 3.0. Average errors for modeled 
versus measured data were 38% for [Chi a] and 37% for ag(400). The high-error values 
were associated with places in the Mississippi River plume where the absorption coefficient 
due to detritus and gelbstoff combined exceeded that for phytoplankton by as much as a 



factor of five. 


c) Data/ Analysis/Interpretation: 
cl): Data collection methods 

The data used to develop and test the algorithm were collected during June 1993 on 
the COLOR cruise at stations shown in Fig. 1. The data taken include spectral remote- 
sensing reflectance (R rs (X)), chlorophyll a and pheopigment a concentrations ([Chi a] and 
[pheo a]), and spectral absorption coefficients for phytoplankton, detritus, and gelbstoff 
(a^(X), a d (X), and a g (X)). 

R r TXt measurements: 

R rs (X) was measured from 390 to 1100 nm using a Spectron Engineering SE-590 
hand-held spectrometer that was modified by adding a 200 /xm slit and a vertical polarizer. 
The polarizer limits the Fresnel reflectance of skylight to values less than 2.1%, which is the 
minimum reflectance value for an unpolarized instrument (Jerlov 1968). This decreases 
uncertainties due to sea-surface waves in the instrumental field of view. An identically 
equipped SE-590 was used simultaneously without a polarizer to evaluate the polarization 
effects on R^X). Spectral ratios of the polarized versus non-polarized data were within 3%, 
indicating that spectral effects due to polarization were not significant. 

Three radiance spectra must be taken by the SE-590 in order to determine R,. S (X). 
These three measurements were taken rapidly in order and the process repeated twice to yield 
three sets of three radiance spectra at each station. First, a water-viewing measurement was 



made at about 25° from nadir, typically at a 90° azimuthal angle from the solar plane. The 
second measurement was of the sky at a 25° zenith angle in the same azimuthal plane as the 
first measurement. The third measurement was of a horizontal Spectralon reflectance 
standard panel with a reflectance of about 10%. It is held by one scientist sighting on the 
horizon across the panel while a second scientist measures the reflected radiance at the same 
orientation as for the water-viewing measurement. The bi-directional reflectance distribution 
of the reflectance standard should be ascertained so that changes in illumination geometry can 
be corrected if necessary. If the panel is Lambertian, xL p (X)/p(X) = E d (X), where Lp is the 
measured panel radiance, p(X) is the irradiance reflectance of the panel, and E d is the 
down welling irradiance. If the panel is not Lambertian, the distribution function 
Q p (0 o ,<f> o ;0,4>) = p(X)E d (X;0 o ,4>o)/L p (X;0,tf>) for the panel should be used instead of x. 

R^X) is the ratio of upwelling water-leaving radiance to the downwelling irradiance 
incident upon the sea surface. It can be calculated from the three radiance measurements via 

* INJjX) - rfX) N s q)] 

N p 

where N w = aL w , N s = aL s , and N p = aL p = opE d /x. N is digital counts, r f is the 
Fresnel reflectance of the sea surface, L is measured radiance, a is the responsivity of the 
SE-590 to radiance, and the subscripts ’w’, ’s’, and ’p’ refer to water, sky, and reflectance 
panel measurements. Note that this approach is not dependent on the instrumental calibration 
factor, a, which divides out. It is only dependent upon the panel calibration factors, which 
remain stable unless the panel surface becomes contaminated. 







pi gment measurements: 

[Chi a] and [pheo a] were determined fluorometrically on a Turner Designs 10-AU- 
005 fluorometer. The method used is described in Holm-Hansen and Reimann (1978), and it 
is an improvement of the method described in Strickland and Parsons (1972). 

a ^(X) and a d (X) measurements; 

The absorption coefficients for particles and detritus, ap(X) and a d (X), were 
determined via the filter pad transmission technique (Mitchell and Kiefer 1988), and the 
methanol extraction technique (Kishino 1985; Roesler et al. 1989). The absorption 
coefficient for phytoplankton, a ^(X), is simply a^X) - a d (X). 

a g (XI measurements: 

The absorption coefficient for gelbstoff, a g (X), was measured at Hg-emission 
wavelengths (365, 404, 436, 546, and 578 nm) using a folded-path 1 m spectrophotometer. 
The regulated Hg arc source was evaluated for each sample measurement, using a calibrated 
reflectance target inserted in the optical path just in front of the sample chamber entrance 
window. The seawater sample was filtered immediately after collection using a 0.2 jim 
Gelman Supor 200 membrane filter, and was refrigerated until the end of the day, when all 
samples collected that day were run consecutively. A deionized water (DIW) blank was run 
before and after the group of samples. DIW and samples were introduced into the 
instrument by pumping through a 0.2 /xm pore size Whatman Polycap 36AS cartridge filter. 
A Spectron Engineering SE-590 spectrometer was used as the sensor. Complete details of 
the a g (X) measurement procedure will be presented at the AGU/ASLO Ocean Sciences 



meeting in February, 1994. 


c2): Al gorithm Development 


R r s model; 

The optical and bio-optical algorithms for SeaWiFS are based on an R,. S (X) model 
which is given by the following general equations (Lee et al. 1992; Lee et al. 1994): 
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R rs is remote-sensing reflectance, b b is the total backscattering coefficient, Q is the upwelling 
irradiance-to-radiance ratio, E^/Ly, a is the total absorption coefficient, and the subscripts 
’w’ and ’p’ refer to water and particles. The lead coefficient 0.176 in Eq. 2 includes effects 
of the air-sea interface, b bw (X) is known (Smith and Baker 1981), and the average value of 
Q W (X) is 3.3 ster' 1 (Lee et al. 1992). X = b bp (400)/Q p (400) and is an indication of the 
magnitude of b bp (X), while Y provides its spectral shape. This Y term is especially 
important near river mouths where extremely fine clays and flocculates can remain after 
heavier particles have settled out of the sediment plume. 



The behavior of X and Y will change with different water types, and we need to 
know a priori a means of classifying water types based simply on season, latitude, and 
R,. S (X). The two Gulf of Mexico data sets include hyper-spectral measurements of R,. S (X) and 
total absorption. X and Y were determined for each station by inverting the model expressed 
in Eqs. 2-4. Expressions were then sought for both X and Y as functions of R,. s at one or 
more of the SeaWiFS wavebands. It was found that X covaries with R rs (670), and linear 
regression on the summer data set yields X = 0.0000328 -I- 3.485 • R rs (670) (n = 15, 
r 2 = . 95 ; see Fig. 2). However, the 670 nm radiance band on SeaWiFS is likely to have too 
low a signal-to-noise ratio to be useful in clearer waters, and we restrict the use of this 
regression equation to waters where R,. s (670) > 0.0008. For clearer waters we use X = 
b bp (400)/Q p (400) and b bp (400) = 0.0036 • [Chi a] 0 25 , which is a function developed using 
the approach for estimating b^ at other wavelengths found in Gordon et al. (1988). [Chi a] 
is estimated via a CZCS-type two-wavelength algorithm determined by linear regression of 
log[Chl a ] vs. log[R rs (443)/R rs (555)] for the same data set, and Q p (400) is set to 4.0 ster' 1 , 
which is an estimate based on near-surface upwelling irradiance and radiance data for the 
Gulf of Mexico provided by W. J. Rhea and C. Davis. Our equations for X for the Gulf of 
Mexico in the summer are summarized below. 

for RJ670) > 0.0008: X = 0.0000328 + 3.485 «„(670) (5) 


for R rs (6 70) < 0.0008: X = 0.0036 [Chi a] 025 /4.0 (6) 



where [Chi a] = 1.71 
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Y was found to covary with the spectral ratio Rj. s (443)/R rs (490) . Note that this ratio 
is not completely independent of absorption by phytoplankton and gelbstoff as both absorb 
more at 443 nm than at 490 nm. For offshore stations this ratio did provide the tightest 
correlation (r 2 > 0.9) with Y of any spectral ratio, and linear regression was used to 
estimate Y for these stations (Fig. 3). For more turbid waters near the river, however, it is 
expected that detritus and fine suspended sediments will produce a slightly blue-rich 
backscatter. We found that a d (412) covaried with R,. s (670), and developed an expression for 
Y as a function of R rs (670). The equations used to calculate Y for the Gulf of Mexico in the 
summer are summarized below. 
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The inverse correlation of Y with [Chi a ] implicit in Eq. 8 suggests that oligotrophic waters 
have smaller particles as has been discussed by others (Gordon and Morel 1983 and 
references cited). Also, the slope value of 3.75 in Eq. 8 was only slightly larger than that 
found for springtime, suggesting that summer backscattering is bluer than in the springtime, 
probably due to the presence of smaller phytoplankton and greater bacterial and detrital 


fractions in the summer. 



Relationships such as these for X and Y are needed on a global basis for a truly 
mature global algorithm to evolve. We believe that sea-surface temperature (SST) can help 
with the parameterization since Y is expected to increase with reductions in nutrients or 
increases in SST. 

Optical algorithm: 

Using spectral ratios of R,. S (X) reduces the uncertainties in transmissivity through the 
air-sea interface that arise due to sea foam and wave-facet angles, 
into Eq. 2, rearranging, and forming a spectral ratio yields 

RJ+l [*U*,)/3.3 + X (400/ r ] 

[^(A.^/3.3 + X (400/ X ; ) y ] 


Substituting Eqs. 3 and 4 
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At this point, the model is still strictly optical in nature (i.e., pigment concentration is not 
involved), and spectral ratios of total absorption coefficients can be determined from 
measurements of R rs (X), given appropriate expressions for X and Y. 

Bio-optical algorithm: 

The bio-optical portion of the algorithm is derived by expanding the total absorption 
coefficient into its constituent parts: 


a(X) = + <*J£) + a g W 
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The subscripts ’w\ ’d’, and ’g’ represent components due to water, phytoplankton, 

detritus, and gelbstoff. If we can develop expressions for a^(X), a d (X), and a g (X) as 
functions only of [Chi a\ and a g (400), then we can combine Eqs. 10 and 11 and use three 
wavelength bands to get a system of two equations in two unknowns. 

The phytoplankton term can be written a^X) = a^(X)[Chl a], where a^(X) is the 
chlorophyll a-specific absorption coefficient for phytoplankton. An empirical relationship 
between nutrient availability (which correlates with chlorophyll concentration) and cell size 
(which affects pigment packaging) seems to hold for much of the ocean (Herbland et al. 

1985; Carder et al. 1986), and its effect on a^(X) has been suggested (Carder et al. 1991) to 
have the form 

= a oW ex P [ a i(^) tank [ a 2 M toCOjCA.) [Chi a])] ] < 12 ) 

The equation parameters a 0 (X) - a 3 (X) can be determined for a given region and season from 
measured a^(X) and [Chi a] data. a^X) is the most important of these parameters, as it is 
directly proportional to a^(X). For simplicity, only a 0 (X) and a,(X) were used to 
parameterize a^(X), with a 2 (X) and a 3 (X) held constant. Fig. 4 shows the measured data and 
the modeled curve for X = 443 nm for the summer cruise, and the parameters for all 
wavelengths are listed in Table 1. The summer a^X) values shown are about twice the 
spring values, most likely due to smaller cells, less [Chi d\ needed per cell due to higher 
light, and more photo-protective pigments required in the summer. 

a d (X) and a g (X) can both be fit to a curve of the form a x (X) = a x (400) exp(-S x (X-400)) 



where the subscript ’x’ refers to either ’d’ or ’g’ (Roesler et al. 1989; Carder et al. 1991). 
Due to this similarity in spectral shape, we eliminate the a^X) term and allow both detrital 
and gelbstoff absorption to be represented by a g (X). For the summer data set, which had 
many stations near the Mississippi River plume, average values of 0.011 nm' 1 and 0.017 
nm' 1 were determined for S d and S g , respectively, and we use an intermediate spectral slope 
value of 0.015 nm" 1 for the combined slope parameter, S. The combined detritus and 
gelbstoff absorption term is thus written 

a (X) = a .(400) e ■° 015 ( l ‘ 40 °) (13) 
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The model equations used in the bio-optical algorithm are derived by inserting Eqs. 

12 and 13 into Eq. 1 1 and substituting the result into both the numerator and denominator of 
Eq. 10. The left hand side of the resulting equation is determined completely by R rs (X) 
inputs, so two of these equations are needed to solve for both [Chi a] and ag(400). Based on 
the shapes of the absorption curves for chlorophyll and gelbstoff, equations using the spectral 
ratios 412/443 and 443/555 should provide the best separation of the two. Signal-to-noise 
considerations may make it necessary to switch to other bands, but the philosophy behind the 
algorithm will remain the same. Using the spectral ratios suggested above yields the 
following equations: 
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These equations are solved via look-up tables (LUTs). First, fifty-element arrays are 
created for both [Chi a] and a g (400). The right hand sides of Eqs. 14 and 15 are then solved 
for each combination of [Chi a ] and a g (400) array values, and the resulting pairs of total 
absorption ratios, a(443)/a(412) and a(555)/a(443), are tabulated in a 2-dimensional array 
(i.e., the LUT). The left hand sides of the equations are evaluated for each input spectrum 
of measured R rs (X), yielding modeled total absorption ratios. The LUT is then searched for 
absorption ratios which are just above and just below the modeled values. Modeled [Chi a ] 
and a g (400) values are then calculated by interpolating between the LUT values of [Chi a] 
and a g (400) that correspond to the LUT absorption ratios determined in the previous step. 
This process can be visualized in Fig. 5 which depicts the LUT as a nomogram of absorption 
ratios, along with the input data points for the summer data set. 



d): Discussion 

The optical portion of the summer algorithm was tested against measured total 
absorption coefficients. Fig. 6a shows the algorithm results versus measured values for 
a(443)/a(412), while those for a(555)/a(443) are shown in Fig. 6b. The algorithm resulted in 
average errors of 7.6% and 11.5%, respectively. These low errors are not surprising, since 
the X and Y parameters used were determined in part from the measured total absorption 
coefficients. 

The bio-optical portion of the summer algorithm was tested against measured values 
of [Chi a] (Fig. 7a) and measured a g (400) (Fig. 7b). This algorithm resulted in average 
errors of 37.9% and 36.6%, respectively. 

The accuracy of estimates of [Chi a] can be limited by many factors, especially at 
high concentrations. Such factors include low signal-to-noise of the measured radiance, high 
absorption by detritus or suspended sediment poorly predicted variations in a^(X), S, and 
high patchiness in the river plumes. 

Errors resulting from low signal-to-noise in L w cannot be easily corrected. Such 
pixels must be flagged and the data processed with an alternate algorithm using wavelengths 
with a stronger signal. Such algorithms will lack the spectral finesse to as accurately 
separate gelbstoff from pigments, so an accuracy flag should be attached to such pixels. 

Errors resulting from excessive detritus or suspended sediments can occur near river 
mouths and other coastal areas. For example, some of the summer cruise stations in the 
Mississippi River plume had a d (412) values larger than 0.4 m' 1 , which was often larger than 
a^(412) and nearly as large as a g (412). This has significant ramifications at 555 nm since 
detrital absorption at that wavelength can be as large or larger than gelbstoff or water 



absorption. In the future, a d (X) and a g (X) should be treated separately in such areas. As a 
preliminary investigation, we found that a d (412) correlated well with 1^(670), and if a d (412) 
can be determined from R rs (670), a d at other wavelengths can be calculated via a d (X) = 
a d (412) exp(-0.011(X-412)). This type of correction will be implemented in the next version 
of the algorithm. Its implementation will, however, require a three-dimensional LUT which 
will increase processing time, so it should only be used in regions with well defined detritus 
problems. 

Unpredictable variations in a^(X) affect [Chi a] retrievals directly since a^(X) is a 
measure of light absorption per unit [Chi a], and these variations can be caused by many 
factors. Revisiting Fig. 4, note the scatter of the measured a^(443) values at the high [Chi 
a] end. The largest deviations of the measured values from the modeled line tend to occur at 
stations where a d (412) > a^(412). Some of the variation is likely also due to non-monotonic 
transitions between dinoflagellate- and diatom-dominated phytoplankton assemblages with 
increasing [Chi a]. This is so because the modeled a^(X) curve is expected to be consistent 
with a smooth, monotonic transition. In addition, light adaptation by phytoplankton can 
effect their pigment packaging and, consequently, a^(X). In fact, this problem appears to be 
the limiting factor in improving algorithm accuracy, at least for these waters in the summer. 


e): Future Actions 

There may be some hope, however, when dealing with changes in a^(X) that are 
associated with nutrient enrichment due to upwelling or vertical mixing. The empirical 
problem is knowing exactly when and how much to vary a^X) and other parameters of Table 



1 to blend the model seamlessly from one season and geographical location to another. 

Species succession and pigment packaging are influenced to a large extent by the nutrients 
available (e.g. no silica means no diatoms), and a satellite-based means of estimating nutrient 
availability would provide a cue for transitioning parameters such as a^X) on a seasonal or 
geographic basis. 

Satellite-derived sea-surface temperature (SST) may provide such a means. Increases 
in SST can indicate increased stability, reduced mixing across the thermocline, and reduced 
nutrients. Decreases in SST can indicate upwelling and/or overturn of the upper layer, 
reduced stability, and increased nutrients. A model relating SST to nutrient (N, P, Si) 
availability on a global basis has been developed using the NODC nutrient, sigma-t and 
temperature data base (Kamykowski and Zentara 1986; Kamykowski 1987). When SST 
exceeds a critical value for a given region of the ocean, a specific nutrient is predicted to be 
limited. Given this model, a species reliant on an influx of the limited nutrients from across 
the nutricline can be expected to be succeeded by phytoplankton ensembles more reliant on 
other nutrients or upon recycled nutrients. Nutrient-availability estimates, not surprisingly, 
would also play a key role in improving the accuracy of primary-production models 
(Kamykowski 1987; Gregg and Walsh 1992). 

A temperature-coupled bio-optical model will identify regions of strong SST 
anomalies and provide a means to make first-order adjustments to the algorithm should that 
be desirable for future versions of the algorithm. Such algorithms require SST measurements 
from AVHRR, which is beyond the scope of our present task for developing algorithms for 
SeaWiFS. We do hope, however, that collaboration with D. Kamykowski will help in the 
initial global parameterization of the model and its seasonal transition from one state to 



another. 
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Table 1. Parameters for the a^(X) vs. [Chi a] relationship (Eq. 12) for data taken on the 
COLOR cruise, June 1993. Only the 412, 443, and 555 nm parameters are used in the 
algorithm described in this paper. 


wavelength 


a i 

a 2 

a 3 

412 

0.040 

0.95 

-0.5 

1.0 

443 

0.060 

0.90 

-0.5 

1.0 

490 

0.039 

0.81 

-0.5 

1.0 

510 

0.026 

0.60 

-0.5 

1.0 

555 

0.008 

0.75 

-0.5 

1.0 
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Fig. 2. X vs. 1^(670) regression (n = 16, ^=0.95). 
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Fig. 4. a^(443) vs. [Chi a]. The dotted line is determined by the parameters in Table 1 and 





Fig. 5. Nomogram depicting the look-up table (LUT) of absorption ratios (dots) and the input 
absorption ratios (diamonds) calculated from R rs data. The dotted lines connect LUT absorption 
ratios that correspond to constant [Chi a] or a g (400). [Chi a\ and a g (400) are calculated from 
the input absorption ratios via bi-linear interpolation of the input points between the four LUT 
grid points that surround it. 
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Fig. 6. Optical algorithm results. Measured vs. modeled absorption ratios for a) 
a(443)/a(412), and b) a(555)/a(443). 
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Fig. 7. Bio-optical algorithm results. 

Measured vs. modeled values for a) [Chi a\, and b) 





